Brandvain and Coop. 

Sperm dependent female meiotic 

drive. 

Model I . Taditional female drive 



Set up 

(★Allele and Genotype frequencies*) 

ClearAll [ "Global^ * " ] 

fA = 1 - fB; 

f AA = f A * 2 + f A f B x ; 

fAB = 2 fA fB (1 - x) ; 

f BB = f B A 2 + f A f B x ; 

Drive 

fAA Drive = FullSimplify [fA ( f AA + fAB (1 - d) ) ] ; 

fAB Drive = FullSimplify [f A (f BB + fAB d) + fB ( f AA + fAB (1 - d) ) ] ; 
fBB Drive = FullSimplify [fB (fAB d + f BB ) ] ; 

Selection 

wAA = 1; wAB = 1 - hs; wBB = 1 - s; (*genotypic fitnesses*) 

W = FullSimplify [fAA Drive wAA + fAB Drive wAB + fBB Drive wBB] ; (*mean fitness*) 
fAA Sel = FullSimplify [ (fAA Drive * wAA) / w] ; 
fAB Sel = FullSimplify [ (f AB Drive * wAB) / w] ; 
fBB Sel = FullSimplify [ (f BB Drive wBB) / w] ; 

fA Se i = FullSimplify [fAAgei + fAB Se i / 2] ; 
fB sei = FullSimplify [fBB Sel + fAB Sel / 2] ; 
AfA = FullSimplify [fA Se i - fA] ; 
AfB = FullSimplify [fB Se i - fB] ; 



Analysis 

For all analyses we assume HWE (i.e. x = 0), we relax this assumption in numerical iterations 
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Invasion 

scaledDeltaDriveWhenRare = (FullSimplify [Af B / (fB * (1 - fB) ) /• x -» 0] /. fB -> 0) 

1 

- (- 1 - 2 d (- 1 + hs) - hs) 
2 

Noting that the change in frequency of a standard driver is independent of s when rare, we solve for the 
value of hs required for invasion. 

invasionCrit = Solve [scaledDeltaDriveWhenRare == 0, hs] 

1 + 2 d 



hs ^ " II 

1 + 2 d J J 



Fixation 

scaledDeltaDriveWhenCommon = (FullSimplify [AfB / (fB * (1-fB)) /. x->0] /. fB -> 1) 

l + hs + 2 (d ( - 1 + hs) - 2hs + s) 



2 (-1 + s) 

We solve for fixation conditions when recessive or not. 

f ixationCritRecessive = Solve [scaledDeltaDriveWhenCommon == 0, s] / . hs -> 0 
f ixationCritNotRecessive = Solve [scaledDeltaDriveWhenCommon == 0, s] 

s -> - (-1 + 2 d) JJ 
s -> - (-l + 2d+3hs-2dhs)}} 

Equilibrium 

We identify the equilibrium frequency of the standard driver 

eqfB = Solve [ (FullSimplify [ (AfB / . x -> 0) ] ) == 0, fB] [ [4] ] 



r 8dhs-4ds + A/-4 (1 - 2 d + hs + 2 d hs) (-4 hs + 8 dhs + 2 s - 4 d s) + (-8dhs + 4ds) 2 , 

f B -> 

1 2 (-4hs + 8dhs + 2 s-4ds) J 

We plot this equilibrium frequency assuming a recessive fitness cost. 
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ContourPlot [If [fB > 1, 1/0, If [fB < 0, 1/0, fB] ] /. eqfB /. hs 0, {d, .5, 1}, 
{s , 0, 1} , PlotLegends -> BarLegend [Automatic, LegendLabel -> " f B* " ] , PlotLabel 

"Equilibrium frequency of a standard driver with a recessive fitness cost", 
FrameLabel -> {"d (drive of the B allele)" , 

"s (selection against driver, assuming recessive fitness cost)"}] 

Equilibrium frequency of a standard driver with a recessive fitness cost 
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dO (drive of the B allele) 

FullSimplify [f AA Drive + fAB Drive + fBB Drive ] 
1 

Model 2. Female drive depends on sperm haplotype 
(single pleitropic locus) 

The B allele is transmited with probability, d, in heterozygous females when fertilized by B-bearing 
sperm. 

x represents the deviation from Hardy - Weinberg Equilibrium 



Setup 

(★Allele and Genotype frequencies*) 

ClearAll [ "Global^ * " ] 

f A = 1 - fB; 

f AA = f A A 2 + f A f B x ; 

fAB = 2 fA fB (1 - x) ; 

f BB = fB A 2 + fAfBx; 
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Drive 

(★Genotype frequencies after drive*) 
fAA Drive = FullSimplify [fA (fAA + fAB / 2) ] ; 

fAB Drive = FullSimplify [fB (fAA + fAB * (1 - d) ) + f A (fAB / 2 + f BB) ] ; 
fBB Drive = FullSimplify [fB (fAB d + f BB ) ] ; 

Selection 

wAA = 1; wAB = 1 - hs; wBB = 1 - s; (*genotypic fitnesses*) 

W = FullSimplify [fAA Drive wAA + fAB Drive wAB + fBB Drive wBB] ; (*mean fitness*) 

fAA Sel = FullSimplify [ (f AA Drive * wAA) / w] ; 

fAB Sel = FullSimplify [ (f AB Drive * wAB) / w] ; 

fBB Sel = FullSimplify [ (fBB Drive wBB) / w] ; 

fA sei = FullSimplify [fAA Sel + fAB Sel / 2] ; 

fB Se i = FullSimplify [fBBgei + fAB Se i / 2] ; 

AfA = FullSimplify [fA Sel - fA] ; 

AfB = FullSimplify [fB Se i - fB] ; 



Analysis 

Note, we assume no deviation from Hardy-Weinberg [i.e. x=0] for all analytical results, and therefore 
these answers are approximations. In the supplamentary material we show thats results of exact recur- 
sions are remarkably consistant from these approximate analystical solutions. 

Assuming the cost of drive is fully recessive [i.e. hs is zero] 
Invasion 

AfBinvade = (FullSimplify [AfB / .hs->0/.x->0]/fB A 2/. fB 0) 

1 

- (-1 + d (2 - 4 s) ) 
2 

spermDepReceesivelnvade = Solve [AfBinvade == 0, s] 

{{-—}} 

plotInvasion4spermDepRecessive = 

Plot[s /. spermDepReceesivelnvade [[1]]/ {d, .5, 1} , PlotStyle -> {Black, Thick}]; 

plotRelChange4RarespermDepRecessive = ContourPlot [ {AfBinvade} , {d, 0.5, 1}, 
{s, 0, 1}, PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AB" ] , 
FrameLabel -> {"d (drive)", "s (selection)"}, 

PlotLabel -> "Invasion of recessive sperm dependent driver"]; 
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Show [plotRelChange4RarespermDepRecessive, plotInvasion4spermDepRecessive] 

Invasion of recessive sperm dependent driver 

scaled AB 

I OA 
0 

- -0.4 

r -0.8 

I 



0.5 0.6 0.7 0.8 0.9 1.0 

d (drive) 

Fixation 

Af Bf ix = FullSimplif y [FullSimplif y [Af B / . hs -» 0 / . x -» 0] / f A] / . f B 1 

-l+2d-2s 
2 - 2 s 

spermDepReceesiveFix = Solve [AfBfix == 0, s] 

{H 2 ( - 1+2d) }} 

(s / . spermDepReceesiveFix [ [1] ] ) 

1 

- (-l + 2d) 
2 

plotFixation4spermDepRecessive = 

Plot[s /. spermDepReceesiveFix [[1]], {d, .5, 1} , PlotStyle -> {Red, Thick}]; 

(*Note we artificially rescaled z to be -.1 for all negative values*) 
plotRelChange4CommonSpermDepRecessive = 

ContourPlot [If [s > (s/. spermDepReceesiveFix [[1]])/ - .1/ AfBfix] , {d, 0.5, 1} , 
{s, 0, 1} , PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AB" ] , 
FrameLabel ^ {"d (dirve)", "s (selection) " } , 

PlotLabel -> "Fixation of recessive sperm dependent driver"]; 




6 | MathematicaAppendix.nb 



Show [plotRelChange4CommonSpermDepRecessive, plotFixation4spermDepRecessive] 

Fixation of recessive sperm dependent driver 
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d (dirve) 



Bistability Point 

FBbistabSpermDepReceesive = Solve [FullSimplify [AfB /.hs ->()/• x->0] == 0 , fB] [ [4] ] 

c l-2d+4ds, 



bistab = ContourPlot [ (If [fB <0, 0, If[fB>l, 1, fB] ] ) /. FBbistabSpermDepReceesive, 
{d, .5 r 1} r {s , 0, 1} , PlotLegends -> BarLegend [Automatic, LegendLabel -> 11 fB*" ], 
FrameLabel -> {"d (dirve)", "s (selection)"}, PlotLabel -> 
"Threshold frequency for fixation of recessive self -promoting driver"]; 



MathematicaAppendix.nb \ 7 



Show [bistab, plotFixation4spermDepRecessive, plotInvasion4spermDepRecessive] 

Threshold frequency for fixation of recessive self-promoting driver 

fB* 

0.9 
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-0.5 
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I 
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d (dirve) 

Assuming the cost of drive is not fully recessive [i.e. hs is 
nonzero] 

Invasion 

Note with any heterozgous cost (i.e. hs > 0) a self - promoting driver cannot invade 

FullSimplif y [FullSimplif y [Af B /. x -> 0] / fB] / • fB -> 0 

-hs 

Fixation 

AfBfix = FullSimplify [FullSimplif y [FullSimplif y[AfB / . x -> 0] / fA] / • fB -> 1] 

I + 2d (-1 + hs) - 3 hs + 2 s 

2 (-1 + s) 

spermDepNotReceesiveFix = Solve [AfBfix == 0, s] 

jjs -> - (-l + 2d+3hs-2dhs)JJ 

spermDepAddFix = Solve [AfBfix == 0 / . hs -» s / 2, s] 

ff 2 (-l + 2d) n 

S ^ 

II 1 + 2 d J J 




I MathematicaAppendix.nb 



plotspermDepAddFix = 

Plot[s /. spermDepAddReceesiveFix , {&, .5, 1} , PlotStyle -> {Red, Thick}]; 

Bistability Point 

FBbistabSpermDepNotReceesive = Solve [FullSimplify [AfB /. x -> 0] == 0, fB] [ [3] ] 

(f B -> -l + 2d + 3hs + 2dhs-4ds- 



An Example of a non - recessive driver [Assuming additivity] 

bistab = ContourPlot [ 

(If[fB<0, 0, If[fB>l, 1, fB] ] ) /. FBbistabSpermDepNotReceesive/. hs (s / 2 
{d, .5, 1} , {s , 0, 1} , PlotLegends -» BarLegend [Automatic , LegendLabel -> " f B* " ] 
FrameLabel -> {"d (dirve)", "s (selection) " } , PlotLabel -> 
"Threshold frequency for fixation of recessive self -promoting driver"]; 

Show [bistab, plotspermDepAddFix] 

Threshold frequency for fixation of recessive self -promoting driver 



8hs (-2hs + 4dhs + 2s-4ds) + (1 - 2 d- 3 hs - 2 dhs + 4 d s) 2 



/ 



(2 (-2 hs + 4 dhs + 2 s - 4 d s) ) 



1.0 



0.8 
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d (dirve) 
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Threshold frequency for fixation of recessive self -promoting driver 

1.0 P ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' ^ 



0.8 

Cannot invade 




0.5 0.6 0.7 0.8 0.9 1.0 

d (dirve) 

Threshold frequency for fixation of recessive self -promoting driver 



0.8 

Cannot invade 
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d (dirve) 
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Model 3. Female drive depends on male genotype 
(single pleitropic locus) 

The B allele is transmited with probability, d and dh, in heterozygous females when fertilized BB and AB 
males, respectively. 

x represents the deviation from Hardy - Weinberg Equilibrium 



Setup 

(★Allele and Genotype frequencies*) 

ClearAll [ "Global^ * " ] 

fA = 1 - fB; 

f AA = f A A 2 + f A f B x ; 

fAB = 2 fA fB (1 - x) ; 

f BB = fB A 2 + fAfBx; 

Drive 

(★Genotype frequencies after drive*) 

fAA Drive = FullSimplif y [ f AA ( f AA + fAB / 2) + fAB ( f AA / 2 + fAB (1 - d het ) / 2) ] ; 
fAB Drive = FullSimplif y [ 

f AA (fAB / 2 + f BB) + fAB (f AA / 2 + fAB / 2 + f BB (1 - d hom ) ) + f BB (f AA + fAB / 2) ] ; 
fBB Drive = FullSimplif y [fAB (fAB d het / 2 + f BB d hom ) + f BB (fAB / 2 + f BB) ] ; 

Selection 

wAA = 1; wAB = 1 - hs; wBB = 1 - s; (*genotypic fitnesses*) 

W = Full Simplify [fAA Drive wAA + fAB Drive wAB + fBB Drive wBB] ; (*mean fitness*) 

fAA Sel = FullSimplif y [ (f AA Drive * wAA) / w] ; 

fAB Sel = FullSimplif y [ (f AB Drive * wAB) / w] ; 

fBB Sel = FullSimplif y [ (fBB Drive wBB) / w] ; 

fA sei = FullSimplif y [fAA Sel + fAB Sel / 2] ; 

fB Se i = FullSimplif y [fBB Se i + fAB Se i / 2] ; 

AfA = FullSimplif y [fA Sel - fA] ; 

AfB = FullSimplif y [fB Se i - fB] ; 
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FullSimplify [f AA Drive + fAB Drive + fBB Drive ] 
1 



Analysis 

Analytical example - recessive fitness cost 



Invasion 

invasion4maleDepRecessive = 
Solve [ ( (FullSimplify [ (AfB /. x -> 0 /. hs -> 0) ] / fB"2) /. fB -> 0) ==0, s] 

r r -1 + 2 d het ! ! 



plotiInvasion4maleDepRecessive = 

Plot[s /. invasion4maleDepRecessive /. d het -> d hom , 
{dhom/ «5, 1} , PlotStyle -> {Black, Thick}]; 



Fixation 

f ixation4maleDepRecessive = 
Solve [ (FullSimplify [ (FullSimplify [ (AfB / . x -> 0 / • hs -> 0) / f A] /. fB -> 1) ] ) == 0, s] 

{{s^ X - (-l + 2d hom )}} 

plotFixation4maleDepRecessive = Plot [s /. f ixation4maleDepRecessive /. d het -> d hom , 
{dhom, .5, 1}, PlotStyle -> {Red, Thick}]; 

Show[Plot[0, {d hom , 0.5, 1}, AxesLabel -» {"d (drive) assuming d het =d hom =d" , 

"s (selection against drive homozygotes) 11 } , PlotRange -> {{.5, 1}, {0, .5}}, 
PlotLabel -> " Invasion and fixation conditions for a male dependent driver 11 ], 
plotFixation4maleDepRecessive, plotiInvasion4maleDepRecessive] 

Invasion and fixation conditions for a male dependent driver 

s (selection against drive homozygotes) 
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Invasion and fixation conditions for a male dependent driver 

s (selection) 




Equlibrium 

FBeqMaleepReceesive = 
Solve [FullSimplify [AfB / . x -> 0 / . hs -> 0 / . d het -> d hom ] == 0, fB] [ [4] ] 

r 2(1-2 d hom + 2 s d hom ) , 

1 (-1 + 2 s) (-l + 2d hom ) J 



eqf ig = ContourPlot [ (If[fB<0, 0, If[fB>l, 1, fB] ] ) /. FBeqMaleepReceesive , 

{dhom/ 0.5, 1} , {s , 0, .5}, PlotLegends -> BarLegend [Automatic , LegendLabel -> "fB*"], 
FrameLabel -> { "d (drive) assuming d het =d hom =d" , 

"s (selection, assuming recessive fitness cost)"}]; 

Show [eqf ig, ploti!nvasion4maleDepRecessive, plotFixation4maleDepRecessive] 




0.5 0.6 0.7 0.8 0.9 1.0 



d (drive) assuming d het =d hom =d 
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0.5 0.6 0.7 0.8 0.9 1.0 



d (drive) assuming dhet=dhom=d 



Model 4. Female drive depends on sperm haplotype 
(two tightly linked, loci in coupling phase) 

We have one locus with two alleles, A (non-driving) and B (traditional driver), as well as a tightly linked 
locus where one allele modifies drive. Assuming no recombination this functions as a third allele, C. 
This tighly-linked locus is on the B background. 

When C increases drive in heterozgous females it fertilizes, it is a drive enhancer (the B+ allele in our 
ms). 

When C decreases drive in heterozgous females it fertilizes, it is a drive suppressor (the B- allele in our 
ms). 
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Setup 

ClearAll [ "Global^ * " ] 

fA = . 

f AA = . 

fAB =. 

fAC =. 

f BB =. 

fBC =. 

fCC =. 

minormod = {dl->dO + e} 

(★assuming the sperm acting modifier additively increases drive by epsilon*) ; 
SUMTOONE = { fA -> 1 - (fB + fC) } ; 
HWE = 

{ f AA -> fA"2, fAB -> 2 f A fB, fAC -> 2 f A fC, f BB -> fB"2, fBC -> 2 f B fC, fCC -> fC"2} ; 
GENOFREQS = { f A -» f AA + fAB / 2 + fAC / 2 , 

f B f BB + fAB / 2 + fBC / 2 , f C f CC + fBC / 2 + fBC / 2 } ; 



Drive 

(*Here we caculate all genotypes after drive. For book- 
keeping purposes we distinguish between reciprocal homozygotes, 
but remove this distinction belowsum them below*) 

AAn = 

FullSimplif y [ f AA * f AA * 1 + f AA * fAB * 1 / 2 + f AA * fAC * 1 / 2 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * (1 - dO) + fAB * fAB * (1 - dO) / 2 + fAB * fAC * (1 - dO) / 2 + 
fAB * f BB * 0 + fAB * fBC * 0 + fAB * f CC * 0 + fAC * f AA * ( 1 - dO) + fAC * fAB * ( 1 - dO) / 2 + 
fAC * fAC * (1 - dO) / 2 + fAC * f BB * 0 + fAC * fBC * 0 + fAC * fCC * 0 + f BB * f AA * 0 + 
f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + 
fBC * fAB * 0 + fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

ABn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 1 / 2 + f AA * fAC * 0 + f AA * f BB * 1 + 
f AA * fBC * 1 / 2 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * ( 1 - dO) / 2 + 
fAB * fAC * 0 + fAB * f BB * (1 - dO) + fAB * fBC * (1 - dO) / 2 + fAB * fCC * 0 + 
fAC * f AA * 0 + fAC * fAB * (1 - dO) / 2 + fAC * fAC * 0 + fAC * f BB * (1 - dO) + 
fAC * fBC * (1 - dO) / 2 + fAC * fCC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 0 + 
fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

ACn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 1 / 2 + f AA * f BB * 0 + 

f AA * fBC * 1 / 2 + f AA * fCC * 1 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (1 - dl) / 2 + 
fAB * f BB * 0 + fAB * fBC * (1 - dl) / 2 + fAB * fCC * (1 - dl) + fAC * f AA * 0 + fAC * fAB * 0 + 
fAC * fAC * (1 - dl) / 2 + fAC * f BB * 0 + fAC * fBC * (1 - dl) / 2 + fAC * fCC * (1 - dl) + 
f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + 
fBC * f AA * 0 + fBC * fAB * 0 + fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + 

fee * f aa * o + fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BAn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * dO + fAB * fAB * dO / 2 + fAB * fAC * dO / 2 + fAB * f BB * 0 + 
fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + fAC * f BB * 0 + 
fAC * fBC * 0 + fAC * fCC * 0 + f BB * f AA * 1 + f BB * fAB * 1 / 2 + f BB * fAC * 1 / 2 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 1 / 2 + fBC * fAB * 1 / 4 + 
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f BC * f AC * 1 / 4 + f BC * f BB * 0 + f BC * f BC * 0 + f BC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BBn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * dO / 2 + fAB * fAC * 0 + fAB * f BB * dO + 
fAB * fBC * dO / 2 + fAB * f CC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + 
fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 1 / 2 + f BB * fAC * 0 + 
f BB * f BB * 1 + f BB * fBC * 1 / 2 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 1 / 4 + 
fBC * fAC * 0 + fBC * f BB * 1 / 2 + fBC * fBC * 1 / 4 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BCn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (dl) / 2 + 

fAB * f BB * 0 + fAB * fBC * (dl) / 2 + fAB * fCC * dl + fAC * f AA * 0 + fAC * fAB * 0 + 

fAC * fAC * 0 + fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + 

f BB * fAC * 1 / 2 + f BB * f BB * 0 + f BB * fBC * 1 / 2 + f BB * f CC * 1 + fBC * f AA * 0 + 

fBC * fAB * 0 + fBC * fAC * 1 / 4 + fBC * f BB * 0 + fBC * fBC * 1 / 4 + fBC * f CC * 1 / 2 + 

fee * f aa * o + fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * f bc * o + fee * fee * o + 0] ; 

CAn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + 

fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * dO + fAC * fAB * dO / 2 + fAC * fAC * dO / 2 + 

fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 

f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 1 / 2 + fBC * fAB * 1 / 4 + 

fBC * fAC * 1 / 4 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 1 + 

fCC * fAB * 1 / 2 + fCC * fAC * 1 / 2 + fCC * f BB * 0 + fCC * fBC * 0 + fCC * fCC * 0 + 0] ; 

CBn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + fAB * fBC * 0 + 
fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * dO / 2 + fAC * fAC * 0 + fAC * f BB * dO + 
fAC * fBC * dO / 2 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 1 / 4 + 
fBC * fAC * 0 + fBC * f BB * 1 / 2 + fBC * fBC * 1 / 4 + fBC * f CC * 0 + f CC * f AA * 0 + 
fCC * fAB * 1 / 2 + fCC * fAC * 0 + fCC * f BB * 1 + fCC * fBC * 1 / 2 + fee * fCC * 0 + 0] ; 

CCn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + 

fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * dl / 2 + 

fAC * f BB * 0 + fAC * fBC * dl / 2 + fAC * f CC * dl + f BB * f AA * 0 + f BB * fAB * 0 + 

f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 0 + 

fBC * fAC * 1 / 4 + fBC * f BB * 0 + fBC * fBC * 1 / 4 + fBC * f CC * 1 / 2 + f CC * f AA * 0 + 

fCC * fAB * 0 + fCC * fAC * 1 / 2 + fCC * f BB * 0 + fCC * fBC * 1 / 2 + fCC * fCC * 1 + 0] ; 

(★Genotype frequencies after drive*) 

fAA Drive = FullSimplif y [AAn] ; 

f AB Drive = FullSimplify [ABn + BAn] ; 

fAC Drive = FullSimplify [ACn + CAn] ; 

fBB Drive = FullSimplify [BBn] ; 

f BC Drive = FullSimplify [BCn + CBn] ; 

fCC Drive = FullSimplify [CCn] ; 

(*check, do allele freqs sum to one?*) 

FullSimplify [ 

FullSimplify [fAA Drive + fAB Drive + fAC Drive + fBB Drive + fBC Drive + fCC Drive ] /. HWE /. SUMTOONE] 

1 



16 | Mathematics Appendix, nb 



Selection 



wAA = 1 ; wAC = wAB = 1 - hs ; wBB = wBC = wCC = 1 - s ; 
W= FullSimplify [ 

(wAA f AA Drive + wAB f AB Drive + wAC f AC Dr i ve + wBB f BB Drive + wBC f BC Dr i ve + wCC fCC Dr i ve ) ] ; 
(★Because the C allele arises on the B background we assume 
it has the same impact on individual fitness*) 
FullSimplify [w / . HWE / . SUMTOONE / . hs -> 0] 

l+(fB + fC) (-2 dl fC + 2 dO fB (-1 + fB + fC) - (fB + fC) (f B + f C - 2 dl f C) ) s 

fAA Sel = fAA Drive wAA/W; 
f AB Sel = f AB Drive wAB / W; 
f AC Se i = f AC Drive wAC / W; 
f BB Se i = f BB Drive wBB / W; 
f BC Se i = f BC Drive wBC / W; 
f CC Se i = f CC Drive wCC / W; 

fA Se i = FullSimplify [fAA Se i + (fAB Se i + fAC Se i) / 2] ; 
fB Se i = FullSimplify [fBB Sel + (fAB Sel + fBC Se i) / 2] ; 
fC Se i = FullSimplify [fCCgei + (fAC Se i + fBC Se i) / 2] ; 
AfA = FullSimplify [fA Sel - fA] ; 
AfB = FullSimplify [fB Se i - fB] ; 
AfC = FullSimplify [fC Se i - fC] ; 

(★Check: do genotype freqs after selection sum to one?*) 

FullSimplify [fA Se i + fB Se i + fC Se i] 

1 



Note, we assume no deviation from Hardy - Weinberg for all analytical results, and therefore these 
answers are approximations. In the supplamentary material we show thats results of exact recursions 
are remarkably consistant from these approximate analystical solutions. 

Invasion of standard driver - note the driver always invades when it has a recessive fitness cost. 
Reassuringly, this reproduces model 1. 

invasionStandardDriver = Solve [ 

(FullSimplify [ (AfB /. GENOFREQS /. HWE /. SUMTOONE / . f C -> 0) / f B] / . f B -> 0) == 0, hs] 



Analysis 



Analysis - a standard driver [i.e. C is absent] 




- 1 + 2 dO 



Fixation of standard driver 
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f ixationStandardDriver = Solve [ 

(FullSimplify [ (AfB / f A / . GENOFREQS / . HWE / . SUMTOONE / . fC 0) ] / . fB -» 1) == 0, s 

s -> - (-l + 2d0 + 3hs-2d0hs)}} 



(★fixation of a standard recessive driver*) 
f ixationStandardDriver /. hs -> 0 



s -> - (-1 + 2 dO) JJ 



Equilibrium 

(★Equilibrium frequency of a standard driver]*) 

eqfB = Solve [ (AfB /. GENOFREQS /. HWE /. SUMTOONE /. fC -> 0) == 0, fB] [ [4] ] 

f B -> 8 dO hs - 4 dO s + 



-4 (1 - 2 dO + hs + 2 dO hs) ( - 4 hs + 8 dO hs + 2 s - 4 dO s ) + ( - 8 dO hs + 4 dO s ) z y 
(2 (-4hs + 8d0hs + 2s-4d0s))J 

(★Plot of equilibrium frequency of standard driver assuming full recessivity*) 

ContourPlot [If [fB > 1, 1/0, If [fB < 0, 1/0, fB] ] /. eqfB /. hs 0, {dO, .5, 1}, 
{s , 0, 1} , PlotLegends -> BarLegend [Automatic, LegendLabel -> " f B* " ] , PlotLabel -> 

"Equilibrium frequency of a standard driver with a recessive fitness cost", 
FrameLabel -> {"dO (drive of the B allele)" , 

"s (selection against driver, assuming recessive fitness cost)"}] 

Eq uilibrium frequency of a standard driver with a recessive fitness c ost 

i.o r ^ ^ ' ^ ^ ^ 



fB* 




-0.9 



0.7 



0.5 



0.7 0.8 
dO (drive of the B allele) 



Invasion of sperm acting drive modifier tighly linked with the 
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driver, and on the driving background 

(★change in frequency of the drive modifier when rare and when alleles 
at the drive locus are in drive-viability equilibrium, mutliplied by 
Wbar/fC [this value is always positive and will not influence the sign]*) 
wbarDeltaSpermDrive = FullSimplify [ 

FullSimplify [WAfC / (fC) /. HWE /. SUMTOONE] /. f C -> 0 / . eqfB /. minormod] 
1 

(hs - s) 

2 (1 - 2 dO) 2 (2 hs - s) 

( - 1 - hs - \[Y V (1 + hs - 2 dO (2 + dO (-2 + s) ) ) (2 hs - s) +2d0 (2 ( - 1 + dO) (-1 + hs) + s) ) e 

(★Plotting change this change in frequency when the sperm acting locus is rare, 
increases drive, and when the fitness cost of drive is fully 
recessive. NOTE: This sperm enhancer of drive cannot invade*) 

ContourPlot [ 

If [fB > 1, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDrive] ] / . eqfB / .e->0.01/.hs->0, 
{dO, 0.5, 1}, {s, 0, 1}, 

PlotLegends -» BarLegend [Automatic, LegendLabel -> "scaled AB + " ] , 
PlotLabel -» "A sperm-acting drive enhancer (the B + haplotype) 

cannot invade\n (assuming a recessive fitness cost of drive)", 
FrameLabel -> {"dO (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 

A sperm-acting drive enhancer (the B + haplotype) cannot invade 



(assuming a recessive fitness cost of drive) 




0.5 0.6 0.7 0.8 0.9 1.0 



dO (drive of the B allele) 

(★Plotting change this change in frequency when the sperm acting 
locus is rare and decreases drive, and when the fitness cost of 
drive is fully recessive. NOTE: This sperm suppresor always invades*) 
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ContourPlot [ 

If [fB > 1, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDrive] ] / . eqfB /.e-»-0.01/.hs-»0, 
{d0, 0.5, 1}, {s, 0, 1}, 

PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AB"" ] , 
PlotLabel -> "A sperm -acting drive suppressor (the B" haplotype) 

always invades\n (assuming a recessive fitness cost of drive)", 
FrameLabel -> {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 

A sperm-acting drive suppressor (the B~ haplotype) always invades 



(assuming a recessive fitness cost of drive) 




0.5 0.6 0.7 0.8 0.9 1.0 



dO (drive of the B allele) 
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(★Plotting change this change in frequency when hte extent 
of drive is mild [dO = .6] and the sperm acting locus is rare 
and increase drive. NOTE: This sperm enhancer never invades*) 
ContourPlot [ (If [fB > 1, 1/0, If [ fB < 0 . 0000001 , 0, wbarDeltaSpermDrive] ] /. eqf B / . 
6 -> 0.01 / . dO -» .6) , {s, 0, 1}, {hs, 0, 1}, 
PlotLegends -» BarLegend [Automatic , LegendLabel -> "scaled AB + " ] , PlotLabel -» 
"A sperm-acting drive enhancer (the B + haplotype) cannot invade\n (assuming 
drive-selection equilibrium of a modest driver, dO = .6)" , 
FrameLabel -> { "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 

A sperm-acting drive enhancer (the B + haplotype) cannot invade 
(assuming drive -selection equilibrium of a modest driver, dO = .6) 

uoF J 



scaled AB + 



-0.0001 
-0.0002 
-0.0003 
-0.0004 
- -0.0005 



0.0 ^ 

0.0 




0.2 0.4 0.6 

s (selection against drive homozygotes) 



(★Plotting change this change in frequency when hte extent of 
drive is mild [dO = .98] and the sperm acting locus is rare and 
increase drive. NOTE: This sperm enhancer never invades*) ContourPlot [ 
(If[fB>l, 1/0, If [fB <= 0.001, -0.000000001, wbarDeltaSpermDrive] ] / . eqf B / . 

e -> 0.01 /. dO .98) , {s, 0, 1}, {hs, 0, 1}, 
PlotLegends -» BarLegend [Automatic, LegendLabel -> "scaled AB + " ] , PlotLabel -» 
11 A sperm-acting drive enhancer (the B + haplotype) cannot invade\n (assuming 
drive-selection equilibrium of a strong driver, dO = .98)", 
FrameLabel -> { "s (selection against drive homozygotes)", 
"hs (selection against drive heterozygotes)"}] 



A sperm-acting drive enhancer (the B + haplotype) cannot invade 
(assuming drive -selection equilibrium of a strong driver, dO = .98) 



1.0 



scaled AB + 
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0.0 ^ 

0.0 



0.2 0.4 0.6 

s (selection against drive homozygotes) 




-0.0005 



-0.001 



A sperm-acting drive enhancer (the B + haplotype) cannot invade 
(assuming drive -selection equilibrium of a strong driver, dO = .98) 



8 0.6 



0.4 - 



0.2 - 




0.0 



0.2 0.4 0.6 

s (selection against drive homozygotes) 



0.8 



1.0 



scaled AB + 



-0.0015 



-0.002 
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Replacement of traditional driver by sperm acting drive 
suppressor tighly linked with the driver, and on the driving 
background 

Equilibrium frequency of drive suppresor 

eqfC = Solve [ (FullSimplify [AfC /. HWE /. SUMTOONE /. minormod /. fB 0] ) ==0, fC] [ [4] ] ; 
(★Equilibrium frequency of a linked, coupled, 

drive suppressor when the fitness costs of drive are recessive*) 

ContourPlot [If [fC < 1 && fC > 0, (fC) , 1/0] /. eqfC / .hs-»0/.e-» - .001, 
{dO, .5, 1}, {s, 0, 1}, PlotLegends -> BarLegend [Automatic, LegendLabel -> 11 fB" *" ] , 
PlotLabel -> "Equlibrium freq. of sperm-acting drive suppressor (the B" 

haplotype) \n (assuming an initial complete driver, dO = 1)", FrameLabel -> 
{"d0 (drive of the B allele)", "s (selection against drive homozygotes) 11 } ] 

Equlibrium freq. of sperm-acting drive suppressor (the B~ haplotype) 



(assuming an initial complete driver, dO = 1) 




0.5 0.6 0.7 0.8 0.9 1.0 



dO (drive of the B allele) 

(★Equilibrium frequency of a linked, coupled, 
drive suppressor when drive is complete*) 
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ContourPlot [If [fC < 1 && fC > 0, (fC) , 1/0] /. eqfC / .hs-»0/.e-» - .001 /. dO -> 1, 
{s, 0, 1} , {hs, 0, 1} , PlotLegends -> BarLegend [Automatic, LegendLabel -> " fB" *"] , 
PlotLabel -> "Equlibrium freq. of a sperm-acting drive suppressor (the 

B" haplotype) \n (assuming an initial complete driver, dO = 1)", 
FrameLabel -> { "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 

Equlibrium freq. of a sperm-acting drive suppressor (the B~ haplotype) 
(assuming an initial complete driver, dO = 1) 



fB" * 

-0.9 
0.7 
0.5 
0.3 
0.1 



0.0 0.2 0.4 0.6 0.8 1.0 

s (selection against drive homozygotes) 

Equilibrium frequency of drive/sperm acting suppresor haplotype is often less than that of the 
standard driver it replaces. 



0.8 
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(*Differecce in equilibrium frequency of the B- and B haplotypes*) 
ContourPlot [ 

If [fC < 1 && fB < 1 && fC > 0&& fB > 0, (fC - fB) , 1/0] /. eqfB /. eqfC /. hs -> 0 / . 

e -> - .001, {d0, .5, 1}, {s, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -> 11 fB" * - fB* 11 ] , PlotLabel -> 
"Difference in equlibrium frequency of a drive suppresor (B" haplotype) \n and 
the driver it replaces (assuming a recessive cost of the drive allele)", 
FrameLabel -> {"dO (drive of the B allele)", 
"s (selection against drive homozygotes) " } ] 

Difference in equlibrium frequency of a drive suppresor (B~ haplotype) 
and the driver it replaces (assuming a recessive cost of the drive allele) 



1.0 



fB" * - fB* 



73 




o.o 



0 



-0.0005 



-0.0015 



-0.002 



-0.001 



0.5 



0.6 



0.7 



0.8 



0.9 



1.0 



dO (drive of the B allele) 
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ContourPlot [ 

If [fC < 1 && fB < 1 && fC > 0&& fB > 0, (fC - fB) , 1/0] /. eqfB /. eqfC / . hs -» (s / 4) / . 

e -> - .001, {dO, .5, 1}, {s, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -> 11 fB" * - fB*" ] , PlotLabel -> 
"Difference in equlibrium frequency of a drive suppresor (B" haplotype) \n 
and the driver it replaces (assuming a heterozygote cost, hs = s/4) " , 
FrameLabel -> {"d0 (drive of the B allele)", 
"s (selection against drive homozygotes) " } ] 

Difference in equlibrium frequency of a drive suppresor (B~ haplotype) 
and the driver it replaces (assuming a heterozygote cost, hs = s/4) 



fB" * - fB* 

o 

-0.001 
-0.002 
-0.003 
-0.004 



0.5 0.6 0.7 0.8 0.9 1.0 

dO (drive of the B allele) 

(★change in frequency of a rare traditional driver when the drive- 
sperm suppressor haplotype and the ondriving hplotype 

are at drive selection equilibrium, mutliplied by 
Wbar/fC [this value is always positive and will not influence the sign]*) 
wbarDeltaTradDrive = FullSimplify [ 

FullSimplify [WAfB / (fB) /. HWE /. SUMTOONE] /. f B -» 0 / • eqfC /. minormod] ; 

(★Plotting the change in frequency of a rare traditional 
driving haplotype when the sperm drive suppresor is at drive 
selection balance. ASSUMING the fitness cost of drive is fully 
recessive. NOTE: This sperm enhancer of drive cannot invade*) 
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ContourPlot [ 

If [fC > 1, 1 / 0, If [fC < 0.0001, 0, wbarDeltaTradDrive] ] /. eqfC /. e -± -0.001 /. hs -> 0, 
{d0, 0.5, 1}, {s, 0, 1}, 

PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AB" ] , 

PlotLabel -> "A fixed drive suppresor (B" haplotype) cannot be displaced by \n a 



traditional driver (assuming a recessive fitness cost of drive)", 
FrameLabel -> {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 

A fixed drive suppresor (B~ haplotype) cannot be displaced by 
a traditional driver (assuming a recessive fitness cost of drive) 




1.0 



scaled AB 



-0.000025 



-0.000075 



-0.000125 



- -0.000175 



0.0 



0.5 



0.6 



0.7 



0.8 



0.9 



1.0 



dO (drive of the B allele) 



(★Plotting the change in frequency of a rare traditional 
driving haplotype when the sperm drive suppresor is at drive 
selection balance. ASSUMING a fitness cost of drive in homozygotes 
(hs->s/4) . NOTE: This sperm enhancer of drive cannot invade*) 
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ContourPlot [ 

If [fC > 1, 1/0, If [fC < 0.0001, 0, wbarDeltaTradDrive] ] /. eqfC /. e -0.001 /. 

hs->s/4, {d0, 0.5, 1}, {s, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AB" ] , 
PlotLabel -> "A fixed sperm- acting drive suppresor cannot be displaced 
by a traditional \n (assuming a recessive fitness cost of drive)", 
FrameLabel -> {"d0 (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 

A fixed sperm-acting drive suppresor cannot be displaced by a traditional 
(assuming a recessive fitness cost of drive) 



0.8 - 



0.6 - 



5 °- 4 " 
bo 



0.2 - 




scaled AB 



-0.00002 



-0.00004 



-0.00012 



0.7 0.8 
dO (drive of the B allele) 
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ContourPlot [ 

If [fC > 1, 1 / 0, If [fC < 0.0001, 0, wbarDeltaTradDrive] ] /. eqfC /. e -± -0.001 /. dO -> 1, 
{s , 0, 1} , {hs, 0, 1} , PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AB" ], 
PlotLabel -> "A fixed sperm- acting drive suppresor cannot 

be displaced by a traditional \n (assuming perfect drive)" , 
FrameLabel -> { "s (selection against drive homozygotes) " , 
"hs (selection against drive heterozygotes) " } ] 

A fixed sperm-acting drive suppresor cannot be displaced by a traditional 
(assuming perfect drive) 



scaled AB 

- -0.000025 
-0.000075 

- -0.000125 
-0.000175 



0.0 0.2 0.4 0.6 0.8 1.0 

s (selection against drive homozygotes) 



Model 5. Female drive depends on sperm haplotype 
(two tightly linked, loci in repulsion phase) 

We have one locus with two alleles, A (non-driving) and B (traditional driver), as well as a tightly linked 
locus where one allele modifies drive. This tighly-linked locus is on the A background. Assuming no 
recombination this functions as a third allele, C. 

When C increases drive in heterozgous females it fertilizes, it is a drive enhancer (the A+ allele [not 
discussed in our ms]). 

When C decreases drive in heterozgous females it fertilizes, it is a drive suppressor (the A- allele in our 
ms). 
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Setup 

ClearAll [ "Global^ * " ] 

fA = . 

f AA = . 

fAB =. 

fAC =. 

f BB =. 

fBC =. 

fCC =. 

minormod = {dl->dO + e} 

(★assuming the sperm acting modifier additively increases drive by epsilon*) ; 
SUMTOONE = { fA -> 1 - (fB + fC) } ; 
HWE = 

{ f AA -> fA"2, fAB -> 2 f A fB, fAC -> 2 f A fC, f BB -> fB"2, fBC -> 2 f B fC, fCC -> fC"2} ; 
GENOFREQS = { f A -» f AA + fAB / 2 + fAC / 2 , 

f B -» f BB + fAB / 2 + fBC / 2 , f C f CC + fBC / 2 + fBC / 2 } ; 



Drive 

(*Here we caculate all genotypes after drive. For book- 
keeping purposes we distinguish between reciprocal homozygotes, 
but remove this distinction belowsum them below*) 

AAn = 

FullSimplif y [ f AA * f AA * 1 + f AA * fAB * 1 / 2 + f AA * fAC * 1 / 2 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * (1 - dO) + fAB * fAB * (1 - dO) / 2 + fAB * fAC * (1 - dO) / 2 + 
fAB * f BB * 0 + fAB * fBC * 0 + fAB * f CC * 0 + fAC * f AA * 1 / 2 + fAC * fAB * 1 / 4 + 
fAC * fAC * 1 / 4 + fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + 
f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * 0 + 
fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

ABn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 1 / 2 + f AA * fAC * 0 + f AA * f BB * 1 + 

f AA * fBC * 1 / 2 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * ( 1 - dO) / 2 + fAB * fAC * 0 + 

fAB * f BB * (1 - dO) + fAB * fBC * (1 - dO) / 2 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 1 / 4 + 

fAC * fAC * 0 + fAC * f BB * 1 / 2 + fAC * fBC * 1 / 4 + fAC * f CC * 0 + f BB * f AA * 0 + 

f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + 

fBC * fAB * 0 + fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

ACn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 1 / 2 + f AA * f BB * 0 + 

f AA * fBC * 1 / 2 + f AA * fCC * 1 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (1 - dl) / 2 + 
fAB * f BB * 0 + fAB * fBC * (1 - dl) / 2 + fAB * fCC * (1 - dl) + fAC * f AA * 0 + fAC * fAB * 0 + 
fAC * fAC * 1 / 4 + fAC * f BB * 0 + fAC * fBC * 1 / 4 + fAC * f CC * 1 / 2 + f BB * f AA * 0 + 
f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + f BB * fBC * 0 + f BB * f CC * 0 + fBC * f AA * 0 + 
fBC * fAB * 0 + fBC * fAC * 0 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * f CC * 0 + f CC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BAn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * dO + fAB * fAB * dO / 2 + fAB * fAC * dO / 2 + fAB * f BB * 0 + 
fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + fAC * f BB * 0 + 
fAC * fBC * 0 + fAC * fCC * 0 + f BB * f AA * 1 + f BB * fAB * 1 / 2 + f BB * fAC * 1 / 2 + 
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f BB * f BB * 0 + f BB * f BC * 0 + f BB * f CC * 0 + f BC * f AA * dO + f BC * f AB * dO / 2 + 
f BC * f AC * (dO) / 2 + f BC * f BB * 0 + f BC * f BC * 0 + f BC * fCC * 0 + fCC * f AA * 0 + 
fCC * fAB * 0 + fCC * fAC * 0 + fCC * f BB * 0 + fCC * fBC * 0 + fCC * fCC * 0 + 0] ; 
BBn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * dO / 2 + fAB * fAC * 0 + fAB * f BB * dO + 
fAB * fBC * dO / 2 + fAB * f CC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 0 + 
fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 1 / 2 + f BB * fAC * 0 + 
f BB * f BB * 1 + f BB * fBC * 1 / 2 + f BB * f CC * 0 + fBC * f AA * 0 + fBC * fAB * dO / 2 + 
fBC * fAC * 0 + fBC * f BB * dO + fBC * fBC * (dO) / 2 + fBC * fCC * 0 + fCC * f AA * 0 + 

fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

BCn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * (dl) / 2 + 

fAB * f BB * 0 + fAB * fBC * (dl) / 2 + fAB * fCC * dl + fAC * f AA * 0 + fAC * fAB * 0 + 

fAC * fAC * 0 + fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + 

f BB * fAC * 1 / 2 + f BB * f BB * 0 + f BB * fBC * 1 / 2 + f BB * f CC * 1 + fBC * f AA * 0 + 

fBC * fAB * 0 + fBC * fAC * (dl) / 2 + fBC * f BB * 0 + fBC * fBC * (dl) / 2 + fBC * fCC * dl + 

fee * f aa * o + fee * fAB * o + fee * fAC * o + fee * f bb * o + fee * fBC * o + fee * fee * o + 0] ; 

CAn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + 

f AA * fBC * 0 + f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + 
fAB * fBC * 0 + fAB * fCC * 0 + fAC * f AA * 1 / 2 + fAC * fAB * 1 / 4 + fAC * fAC * 1 / 4 + 
fAC * f BB * 0 + fAC * fBC * 0 + fAC * f CC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * fCC * 0 + fBC * f AA * (1 - dO) + fBC * fAB * (1 - dO) / 2 + 
fBC * fAC * (1 - dO) / 2 + fBC * f BB * 0 + fBC * fBC * 0 + fBC * fCC * 0 + fCC * f AA * 1 + 
fCC * fAB * 1 / 2 + fCC * fAC * 1 / 2 + fCC * f BB * 0 + fCC * fBC * 0 + fCC * fCC * 0 + 0] ; 

CBn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + fAB * fBC * 0 + 
fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 1 / 4 + fAC * fAC * 0 + fAC * f BB * 1 / 2 + 
fAC * fBC * 1 / 4 + fAC * fCC * 0 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + 
f BB * f BB * 0 + f BB * fBC * 0 + f BB * fCC * 0 + fBC * f AA * 0 + fBC * fAB * (1 - dO) / 2 + 
fBC * fAC * 0 + fBC * f BB * (1 - dO) + fBC * fBC * (1 - dO) / 2 + fBC * fCC * 0 + fCC * f AA * 0 + 
fCC * fAB * 1 / 2 + fCC * fAC * 0 + fCC * f BB * 1 + fCC * fBC * 1 / 2 + fCC * fCC * 0 + 0] ; 

CCn = FullSimplif y [ f AA * f AA * 0 + f AA * fAB * 0 + f AA * fAC * 0 + f AA * f BB * 0 + f AA * fBC * 0 + 
f AA * fCC * 0 + fAB * f AA * 0 + fAB * fAB * 0 + fAB * fAC * 0 + fAB * f BB * 0 + fAB * fBC * 0 + 
fAB * fCC * 0 + fAC * f AA * 0 + fAC * fAB * 0 + fAC * fAC * 1 / 4 + fAC * f BB * 0 + fAC * fBC * 1 / 4 + 
fAC * fCC * 1 / 2 + f BB * f AA * 0 + f BB * fAB * 0 + f BB * fAC * 0 + f BB * f BB * 0 + 
f BB * fBC * 0 + f BB * fCC * 0 + fBC * f AA * 0 + fBC * fAB * 0 + fBC * fAC * (1 - dl) / 2 + 
fBC * f BB * 0 + fBC * fBC * (1 - dl) / 2 + fBC * fCC * (1 - dl) + fCC * f AA * 0 + 
fCC * fAB * 0 + fCC * fAC * 1 / 2 + fCC * f BB * 0 + fCC * fBC * 1 / 2 + fCC * fCC * 1 + 0] ; 

(★Genotype frequencies after drive*) 
fAA Drive = FullSimplif y [AAn] ; 
fAB Drive = FullSimplif y [ABn + BAn] ; 
fAC Drive = FullSimplif y [ACn + CAn] ; 
fBB Drive = FullSimplify [BBn] ; 
fBC Drive = FullSimplify [BCn + CBn] ; 
fCC Drive = FullSimplify [CCn] ; 
(*check, do allele freqs sum to one?*) 
FullSimplify [ 

FullSimplify [fAA Drive + fAB Drive + fAC Drive + fBB Drive + fBC Drive + fCC Drive ] / . HWE / . SUMTOONE] 

1 



Selection 
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wAA = wAC = wCC = 1 ; wAB = wBC = 1 - hs ; wBB = 1 - s ; 
W= FullSimplify [ 

(wAA f AA Drive + wAB f AB Drive + wAC f AC Drive + wBB f BB Drive + wBC f BC Drive + wCC fCC Drive ) ] ; 
(★Because the C allele arises on the B background we assume 
it has the same impact on individual fitness*) 
FullSimplify [w / . HWE / . SUMTOONE / . hs -» 0] 

1 + (2 dO (-1 + fB) - fB) fB 2 s 

fAA Sel = fAA Drive wAA/W; 
f AB Sel = f AB Drive wAB / W; 
f AC Se i = f AC Drive wAC / W; 
f BB Se i = f BB Drive wBB / W; 
f BC Se i = f BC Drive wBC / W; 
f CC Se i = f CC Drive wCC / W; 

fA Se i = FullSimplify [fAAgei + (fAB Se i + fAC Se i) / 2] ; 
fB Se i = FullSimplify [fBBgei + (fAB Se i + fBC Se i) / 2] ; 
f C Se i = FullSimplify [fCCsei + (fAC Se i + fBC Se i) / 2] ; 
AfA = FullSimplify [fA Se i - fA] ; 
AfB = FullSimplify [fB Se i - fB] ; 
AfC = FullSimplify [fCgei - fC] ; 

(*Check: do genotype freqs after selection sum to one?*) 

FullSimplify [fA Se i + fB Se i + fC Se i] 

1 

FullSimplify [wfACgei /• HWE /. SUMTOONE] 

2 (-1+ (-1 + dO+dl) fB) fC (-1 + fB + fC) 



Analysis 

Analysis - a standard driver [i.e. C is absent] 

THIS IS THE SAME AS ABOVE, AND WE RENETER IT SIMPLY TO LOAD RESULTS INTO MEM- 
ORY. PLOTS ARE NOT RECREATED 

invasionStandardDriver = Solve [ 

(FullSimplify [ (AfB /. GENOFREQS /. HWE /. SUMTOONE / . fC -» 0) / f B] / . f B -» 0) == 0, 
hs] ; 

f ixationStandardDriver = Solve [ 

(FullSimplify [ (AfB / f A / . GENOFREQS / . HWE / . SUMTOONE / . fC -> 0) ] / . fB -> 1) ==0, 
s]; 

eqfB = Solve [ (AfB / . GENOFREQS / . HWE / . SUMTOONE / . fC -> 0) ==0, fB] [ [4] ] ; 
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Invasion of sperm acting drive modifier tighly linked with the 
driver, on the non-driving background 

(★change in frequency of the drive modifier when rare and when alleles 
at the drive locus are in drive-viability equilibrium, mutliplied by 
Wbar/fC^2 [this value is always positive and will not influence the sign]*) 
wbarDeltaSpermDriveRecessive = FullSimplify [ 

(FullSimplify [w AfC / (fC^2) /. HWE /. SUMTOONE /. eqfB /. hs -> 0] ) /. minormod] 

(-2 dO s + \[Y V (-1 + 2 dO (2 +d0 (-2 + s) ) ) s ) e 
2 (-1 + 2 dO) s 

ContourPlot [ 

If[fB>l, 1/0, If [fB < 0.0001, 0, wbarDeltaSpermDriveRecessive] ] /. eqfB/. 

e->-0.01 / . hs 0, {dO, 0.5, 1}, {s, 0, 1}, 
PlotLegends -» BarLegend [Automatic, LegendLabel -> "scaled AA"" ] , 

PlotLabel -> "A sperm-acting drive suppressor linked\n to the nondriving allele 

(the A" haplotype) \n invades when the cost of drive is fully recessive", 
FrameLabel { "dO (drive of the B allele)", 

"s (selection against driver, assuming recessive fitness cost)"}] 



A sperm-acting drive suppressor linked 
to the nondriving allele (the A~ haplotype) 
invades when the cost of drive is fully recessive 




0.5 0.6 0.7 0.8 0.9 1.0 

dO (drive of the B allele) 



(★change in frequency of the drive modifier when rare and when alleles 
at the drive locus are in drive-viability equilibrium, mutliplied by 
Wbar/fC [this value is always positive and will not influence the sign]*) 
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wbarDeltaSpermDrive = FullSimplify [ 

FullSimplify [WAfC / (fC) /. HWE /. SUMTOONE] /. f C -> 0 / . eqfB /. minormod] 

1 hs (- 1 - hs - 

2 (1 - 2 dO) 2 (2 hs - s) 

\[Y V (1 + hs - 2 dO (2 + dO (-2 + s))) (2hs-s) + 2 dO (2 (-1 + dO) (-l + hs)+s))e 

(★Plotting change this change in frequency when the sperm acting locus is rare, 
increases drive, and when the fitness cost of drive is fully 
recessive. NOTE: This sperm enhancer of drive cannot invade*) 

ContourPlot ^ 

If [fB > 1, 1 / 0, If [fB < 0.0001, 0, wbarDeltaSpermDrive] ] /. eqfB /. e -> -0.01 /. 

hs->s/4, {dO, 0.5, 1}, {s, 0, 1}, 
PlotLegends -> BarLegend [Automatic, LegendLabel -> "scaled AA~ 11 ] , PlotLabel -> 

"Invasion of a sperm-acting drive suppressor linked to the nondriving\n allele 
(the A" haplotype) , assuming the cost of drive is not fully recessive.", 

FrameLabel -> |"d0 (drive of the B allele)", 

"s (selection against driver, assuming het fitness cost, hs= — " Jj 

Invasion of a sperm-acting drive suppressor linked to the nondriving 
allele (the A~ haplotype), assuming the cost of drive is not fully recessive. 




scaled AA 

0.0004 
0.0003 
0.0002 
0.0001 



0.5 0.6 0.7 0.8 0.9 1.0 

dO (drive of the B allele) 
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Model 6. Female drive depends on sperm genotype at 
an unlinked locus 

In addition two our one locus with two alleles, A (non-driving) and B (traditional driver), we have a 
second locus with loci, E and F. We are concerned with the invasion and fixation of the F allele, which 
acts to modify drive when it fertilizes a heterozygote. Drive of the B allele is dO, but changes to d1 when 
fertilized by an 'F' bearing sperm 



Setup 



in[i96]:= f A = . ; f B = . ; f E = . ; f F = . ; 
fAE = fA fE + D; 
fAF = f A f F - D; 
f BE = fB fE - D; 
f BF = fB fF + D; 

randomHaps = { f AAEE -> FullSimplif y [ fAE fAE] , 

fABEE -» FullSimplif y [2 fAE f BE] , fBBEE -> FullSimplif y [ f BE f BE] , 

fAAEF -> FullSimplify [2 fAE fAF] , fABEF -> FullSimplif y [2 (fAE f BF + fAF f BE) ] , 

fBBEF -> FullSimplify [2 f BE f BF] , fAAFF -> FullSimplif y [ fAF fAF] , 

fABFF -> FullSimplify [2 fAF f BF] , fBBFF -> FullSimplif y [ fBF f BF] } ; 

sum2one = {fA ^ 1 - fB, f E -> 1 - fF} 

minormod = {dl->dO + e} 
(★assuming the sperm acting modifier additively increases drive by epsilon*) 

Out[202]= {fA -> 1 - fB, f E -> 1 - fF} 



Drive 

fAAEE Drive = FullSimplify [ 

f AAEE * f AAEE * 1 + f AAEE * fABEE * 1 / 2 + f AAEE * fBBEE * 0 + f AAEE * fAAEF * 1 / 2 + 
f AAEE * fABEF * 1 / 4 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + f AAEE * fABFF * 0 + 
f AAEE * fBBFF * 0 + fABEE * f AAEE * (1 - dO) + fABEE * fABEE * (1 - dO) / 2 + fABEE * fBBEE * 0 + 
fABEE * fAAEF * (1 - dO) / 2 + fABEE * fABEF * (1 - dO) / 4 + fABEE * fBBEF * 0 + 
fABEE * fAAFF * 0 + fABEE * fABFF * 0 + fABEE * fBBFF * 0 + fBBEE * f AAEE * 0 + 
fBBEE * fABEE * 0 + fBBEE * fBBEE * 0 + fBBEE * fAAEF * 0 + fBBEE * fABEF * 0 + 
fBBEE * fBBEF * 0 + fBBEE * fAAFF * 0 + fBBEE * fABFF * 0 + fBBEE * fBBFF * 0 + 
fAAEF * f AAEE * 1 / 2 + fAAEF * fABEE * 1 / 4 + fAAEF * fBBEE * 0 + fAAEF * fAAEF * 1 / 4 + 
fAAEF * fABEF * 1 / 8 + fAAEF * fBBEF * 0 + fAAEF * fAAFF * 0 + fAAEF * fABFF * 0 + 
fAAEF * fBBFF * 0 + fABEF * f AAEE * (1 - dO) / 2 + fABEF * fABEE * (1 - dO) / 4 + 
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f ABEF 


•k 


fBBEE 
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fABEF 
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fAAEF 
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(1 


- dO) / 4 


+ fABEF * fABEF * (1 - dO) / 8 + 


fABEF 
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fBBEF 


k 
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f BBEF 


k 


fABEF 


k 


o 


+ 


fBBEF 


k 


fBBEF 


k 


0 + 


fBBEF * 


fAAFF * 0 


+ 


fBBEF * 
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fABFF 


* 


fBBEE 


k 


0 


+ 


fABFF 


k 


fAAEF 


k 


0 + 
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+ 


fABFF * 
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k 


fBBEE 


k 


0 + 


fBBFF * 


fAAEF k 0 


+ 


fBBFF * 


fABEF * 0 + 


fBBFF 


* 


fBBEF 


k 


0 
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0 + 
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fABFF * 0 


+ 


fBBFF * 


fBBFF * 0 + 0] 



fABEE Drive = FullSimplify [fAAEE * fAAEE * 0 + fAAEE * fABEE k 1 / 2 + fAAEE * fBBEE * 1 + 
fAAEE k fAAEF * 0 + fAAEE k fABEF * 1 / 4 + fAAEE k fBBEF k 1 / 2 + fAAEE k fAAFF * 0 + 
fAAEE k fABFF * 0 + fAAEE k fBBFF * 0 + fABEE k fAAEE * dO + fABEE k fABEE k 1 / 2 + 
fABEE * fBBEE * (1 - dO) + fABEE * fAAEF * dO / 2 + fABEE * fABEF * 1 / 4 + 
fABEE * fBBEF * (1 - dO) / 2 + fABEE * fAAFF * 0 + fABEE * fABFF * 0 + fABEE * fBBFF * 0 + 
fBBEE k fAAEE * 1 + fBBEE k fABEE k 1 / 2 + fBBEE k fBBEE * 0 + fBBEE k fAAEF k 1 / 2 + 
fBBEE k fABEF * 1 / 4 + fBBEE k fBBEF * 0 + fBBEE k fAAFF * 0 + fBBEE k fABFF * 0 + 
fBBEE * fBBFF * 0 + fAAEF * fAAEE * 0 + fAAEF * fABEE * 1 / 4 + fAAEF * fBBEE k 1 / 2 + 
fAAEF * fAAEF * 0 + fAAEF * fABEF * 1 / 8 + fAAEF * fBBEF * 1 / 4 + fAAEF * fAAFF * 0 + 
fAAEF k fABFF * 0 + fAAEF k fBBFF * 0 + fABEF k fAAEE * dO / 2 + fABEF k fABEE * 1 / 4 + 
fABEF * fBBEE * (1 - dO) / 2 + fABEF * fAAEF * dO / 4 + fABEF * fABEF * 1 / 8 + 
fABEF * fBBEF * (1 - dO) / 4 + fABEF * fAAFF * 0 + fABEF * fABFF * 0 + fABEF * fBBFF * 0 + 
fBBEF * fAAEE k 1 / 2 + fBBEF * fABEE * 1 / 4 + fBBEF * fBBEE * 0 + fBBEF * fAAEF * 1 / 4 + 
fBBEF k fABEF * 1 / 8 + fBBEF k fBBEF * 0 + fBBEF k fAAFF * 0 + fBBEF k fABFF * 0 + 
fBBEF k fBBFF * 0 + fAAFF k fAAEE * 0 + fAAFF k fABEE * 0 + fAAFF k fBBEE * 0 + 
fAAFF * fAAEF * 0 + fAAFF * fABEF * 0 + fAAFF * fBBEF * 0 + fAAFF * fAAFF * 0 + 
fAAFF * fABFF * 0 + fAAFF * fBBFF * 0 + fABFF * fAAEE * 0 + fABFF * fABEE * 0 + 
fABFF * fBBEE * 0 + fABFF * fAAEF * 0 + fABFF * fABEF * 0 + fABFF * fBBEF * 0 + 
fABFF * fAAFF * 0 + fABFF * fABFF * 0 + fABFF * fBBFF * 0 + fBBFF * fAAEE * 0 + 
fBBFF * fABEE * 0 + fBBFF * fBBEE * 0 + fBBFF * fAAEF * 0 + fBBFF * fABEF * 0 + 
fBBFF * fBBEF * 0 + fBBFF * fAAFF * 0 + fBBFF * fABFF * 0 + fBBFF * fBBFF * 0 + 0] 

fBBEE Drive = FullSimplify [fAAEE k fAAEE * 0 + fAAEE k fABEE * 0 + fAAEE k fBBEE * 0 + 
fAAEE k fAAEF * 0 + fAAEE k fABEF * 0 + fAAEE k fBBEF * 0 + fAAEE k fAAFF * 0 + 
fAAEE k fABFF * 0 + fAAEE k fBBFF * 0 + fABEE k fAAEE * 0 + fABEE k fABEE * dO / 2 + 
fABEE k fBBEE * dO + fABEE k fAAEF * 0 + fABEE k fABEF k (dO / 4) + 
fABEE k fBBEF * dO / 2 + fABEE k fAAFF * 0 + fABEE k fABFF * 0 + fABEE k fBBFF * 0 + 
fBBEE k fAAEE * 0 + fBBEE k fABEE * 1 / 2 + fBBEE k fBBEE * 1 + fBBEE k fAAEF * 0 + 
fBBEE * fABEF * 1 / 4 + fBBEE * fBBEF * 1 / 2 + fBBEE * fAAFF * 0 + fBBEE * fABFF * 0 + 
fBBEE k fBBFF * 0 + fAAEF k fAAEE * 0 + fAAEF k fABEE * 0 + fAAEF k fBBEE * 0 + 
fAAEF k fAAEF * 0 + fAAEF k fABEF * 0 + fAAEF k fBBEF * 0 + fAAEF k fAAFF * 0 + 
fAAEF k fABFF * 0 + fAAEF k fBBFF * 0 + fABEF k fAAEE * 0 + fABEF k fABEE k (dO / 4) + 
fABEF * fBBEE * dO / 2 + fABEF * fAAEF * 0 + fABEF * fABEF * (dO / 8) + 
fABEF * fBBEF * dO / 4 + fABEF * fAAFF * 0 + fABEF * fABFF * 0 + fABEF * fBBFF * 0 + 
fBBEF k fAAEE * 0 + fBBEF k fABEE * 1 / 4 + fBBEF k fBBEE * 1 / 2 + fBBEF k fAAEF * 0 + 
fBBEF * fABEF * 1 / 8 + fBBEF * fBBEF * 1 / 4 + fBBEF * fAAFF * 0 + fBBEF * fABFF * 0 + 
fBBEF * fBBFF * 0 + fAAFF * fAAEE * 0 + fAAFF * fABEE * 0 + fAAFF * fBBEE * 0 + 
fAAFF * fAAEF * 0 + fAAFF * fABEF * 0 + fAAFF * fBBEF * 0 + fAAFF * fAAFF * 0 + 
fAAFF * fABFF * 0 + fAAFF * fBBFF * 0 + fABFF * fAAEE * 0 + fABFF * fABEE * 0 + 
fABFF * fBBEE * 0 + fABFF * fAAEF * 0 + fABFF * fABEF * 0 + fABFF * fBBEF * 0 + 
fABFF * fAAFF * 0 + fABFF * fABFF * 0 + fABFF * fBBFF * 0 + fBBFF * fAAEE * 0 + 
fBBFF * fABEE * 0 + fBBFF * fBBEE * 0 + fBBFF * fAAEF * 0 + fBBFF * fABEF * 0 + 
fBBFF * fBBEF * 0 + fBBFF * fAAFF * 0 + fBBFF * fABFF * 0 + fBBFF * fBBFF * 0 + 0] 

f AAEF Drive = FullSimplify [fAAEE k fAAEE * 0 + fAAEE k fABEE * 0 + fAAEE k fBBEE * 0 + 
fAAEE * fAAEF * 1 / 2 + fAAEE * fABEF * 1 / 4 + fAAEE * fBBEF * 0 + fAAEE * fAAFF * 1 + 
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f AAEE * f ABFF * 1 / 2 + f AAEE * f BBFF * 0 + f ABEE * f AAEE * 0 + f ABEE * f ABEE * 0 + 

f ABEE * fBBEE * 0 + f ABEE * f AAEF * (1 - dl) / 2 + f ABEE * fABEF * (1 - dl) / 4 + 

fABEE * fBBEF * 0 + fABEE * fAAFF * (1 - dl) + fABEE * f ABFF * (1 - dl) / 2 + 

f ABEE * f BBFF * 0 + fBBEE * f AAEE * 0 + fBBEE * fABEE * 0 + fBBEE * fBBEE * 0 + 

fBBEE * f AAEF * 0 + fBBEE * fABEF * 0 + fBBEE * fBBEF * 0 + fBBEE * fAAFF * 0 + 

fBBEE * f ABFF * 0 + fBBEE * f BBFF * 0 + f AAEF * f AAEE * 1 / 2 + f AAEF * fABEE * 1 / 4 + 

f AAEF * fBBEE * 0 + f AAEF * f AAEF * 1 / 2 + f AAEF * fABEF * 1 / 4 + f AAEF * fBBEF * 0 + 

f AAEF * fAAFF * 1 / 2 + f AAEF * f ABFF * 1 / 4 + f AAEF * f BBFF * 0 + fABEF * f AAEE * (1 - dO) / 2 

fABEF * fABEE * (1 - dO) / 4 + fABEF * fBBEE * 0 + fABEF * f AAEF * ( (1 - dO) / 4 + (1 - dl) / 4) 

fABEF * fABEF * ( (1 - dO) / 8 + (1 - dl) / 8) + fABEF * fBBEF * 0 + fABEF * fAAFF * (1 - dl) / 2 

fABEF * f ABFF * (1 - dl) / 4 + fABEF * f BBFF * 0 + fBBEF * f AAEE * 0 + fBBEF * fABEE * 0 + 

fBBEF * fBBEE * 0 + fBBEF * f AAEF * 0 + fBBEF * fABEF * 0 + fBBEF * fBBEF * 0 + 

fBBEF * fAAFF * 0 + fBBEF * f ABFF * 0 + fBBEF * f BBFF * 0 + fAAFF * f AAEE * 1 + 

fAAFF * fABEE * 1 / 2 + fAAFF * fBBEE * 0 + fAAFF * f AAEF * 1 / 2 + fAAFF * fABEF * 1 / 4 + 

fAAFF * fBBEF * 0 + fAAFF * fAAFF * 0 + fAAFF * f ABFF * 0 + fAAFF * f BBFF * 0 + 

f ABFF * f AAEE * dO + f ABFF * fABEE * (1 - dO) / 2 + f ABFF * fBBEE * 0 + 

f ABFF * f AAEF * (1 - dO) / 2 + f ABFF * fABEF * (1 - dO) / 4 + f ABFF * fBBEF * 0 + 

f ABFF * fAAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 

f BBFF * fABEE * 0 + f BBFF * fBBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * fABEF * 0 + 

f BBFF * fBBEF * 0 + f BBFF * fAAFF * 0 + f BBFF * f ABFF * 0 + f BBFF * f BBFF * 0 + 0] 

f ABEF Drive = FullSimplify [f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * fBBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * fABEF * 1 / 4 + f AAEE * fBBEF * 1 / 2 + f AAEE * fAAFF * 0 + 
f AAEE * f ABFF * 1 / 2 + f AAEE * f BBFF * 1 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * fBBEE * 0 + fABEE * f AAEF * dl / 2 + fABEE * fABEF * 1 / 4 + 
fABEE * fBBEF * (1 - dl) / 2 + fABEE * fAAFF * dl + fABEE * f ABFF * 1 / 2 + 
fABEE * f BBFF * (1 - dl) + fBBEE * f AAEE * 0 + fBBEE * fABEE * 0 + fBBEE * fBBEE * 0 + 
fBBEE * f AAEF * 1 / 2 + fBBEE * fABEF * 1 / 4 + fBBEE * fBBEF * 0 + fBBEE * fAAFF * 1 + 
fBBEE * f ABFF * 1 / 2 + fBBEE * f BBFF * 0 + f AAEF * f AAEE * 0 + f AAEF * fABEE * 1 / 4 + 
f AAEF * fBBEE * 1 / 2 + f AAEF * f AAEF * 0 + f AAEF * fABEF * 1 / 4 + f AAEF * fBBEF * 1 / 2 + 
f AAEF * fAAFF * 0 + f AAEF * f ABFF * 1 / 4 + f AAEF * f BBFF * 1 / 2 + fABEF * f AAEE * dO / 2 + 
fABEF * fABEE * 1 / 4 + fABEF * fBBEE * (1 - dO) / 2 + fABEF * f AAEF * (dO / 4 + dl / 4) + 
fABEF * fABEF * 1 / 4 + fABEF * fBBEF * ( (1 - dO) / 4 + (1 - dl) / 4) + fABEF * fAAFF * dl / 2 + 
fABEF * f ABFF * 1 / 4 + fABEF * f BBFF * (1 - dl) / 2 + fBBEF * f AAEE * 1 / 2 + 
fBBEF * fABEE * 1 / 4 + fBBEF * fBBEE * 0 + fBBEF * f AAEF * 1 / 2 + fBBEF * fABEF * 1 / 4 + 
fBBEF * fBBEF * 0 + fBBEF * fAAFF * 1 / 2 + fBBEF * f ABFF * 1 / 4 + fBBEF * f BBFF * 0 + 
fAAFF * f AAEE * 0 + fAAFF * fABEE * 1 / 2 + fAAFF * fBBEE * 1 + fAAFF * f AAEF * 0 + 
fAAFF * fABEF * 1 / 4 + fAAFF * fBBEF * 1 / 2 + fAAFF * fAAFF * 0 + fAAFF * f ABFF * 0 + 
fAAFF * f BBFF * 0 + f ABFF * f AAEE * (1 - dO) + f ABFF * fABEE * 1 / 2 + f ABFF * fBBEE * dO + 
f ABFF * f AAEF * dO / 2 + f ABFF * fABEF * 1 / 4 + f ABFF * fBBEF * (1 - dO) / 2 + 
f ABFF * fAAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 1 + 
f BBFF * fABEE * 1 / 2 + f BBFF * fBBEE * 0 + f BBFF * f AAEF * 1 / 2 + f BBFF * fABEF * 1 / 4 + 
f BBFF * fBBEF * 0 + f BBFF * fAAFF * 0 + f BBFF * f ABFF * 0 + f BBFF * f BBFF * 0 + 0] 

fBBEF Drive = FullSimplify [f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * fBBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * fABEF * 0 + f AAEE * fBBEF * 0 + f AAEE * fAAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * f BBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * fBBEE * 0 + fABEE * f AAEF * 0 + fABEE * fABEF * dl / 4 + fABEE * fBBEF * dl / 2 + 
fABEE * fAAFF * 0 + fABEE * f ABFF * dl / 2 + fABEE * f BBFF * dl + fBBEE * f AAEE * 0 + 
fBBEE * fABEE * 0 + fBBEE * fBBEE * 0 + fBBEE * f AAEF * 0 + fBBEE * fABEF * 1 / 4 + 
fBBEE * fBBEF * 1 / 2 + fBBEE * fAAFF * 0 + fBBEE * f ABFF * 1 / 2 + fBBEE * f BBFF * 1 + 
f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * fBBEE * 0 + f AAEF * f AAEF * 0 + 
f AAEF * fABEF * 0 + f AAEF * fBBEF * 0 + f AAEF * fAAFF * 0 + f AAEF * f ABFF * 0 + 
f AAEF * f BBFF * 0 + fABEF * f AAEE * 0 + fABEF * fABEE * (dO / 4) + fABEF * fBBEE * (dO / 2) + 
fABEF * f AAEF * 0 + fABEF * fABEF * (dO / 8 + dl / 8) + fABEF * fBBEF * (dO / 4 + dl / 4) + 
fABEF * fAAFF * 0 + fABEF * f ABFF * dl / 4 + fABEF * f BBFF * (dl / 2) + 
fBBEF * f AAEE * 0 + fBBEF * fABEE * 1 / 4 + fBBEF * fBBEE * 1 / 2 + fBBEF * f AAEF * 0 + 
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f BBEF * f ABEF * 1 / 4 + f BBEF * f BBEF * 1 / 2 + f BBEF * f AAFF * 0 + f BBEF * f ABFF * 1 / 4 + 
f BBEF * f BBFF * 1 / 2 + f AAFF * f AAEE * 0 + f AAFF * f ABEE * 0 + f AAFF * f BBEE * 0 + 
f AAFF * f AAEF * 0 + f AAFF * f ABEF * 0 + f AAFF * f BBEF * 0 + f AAFF * f AAFF * 0 + 
f AAFF * f ABFF * 0 + f AAFF * f BBFF * 0 + f ABFF * f AAEE * 0 + f ABFF * f ABEE * dO / 2 + 
f ABFF * f BBEE * (1 - dO) + f ABFF * f AAEF * 0 + f ABFF * f ABEF * dO / 4 + f ABFF * f BBEF * dO / 2 + 
f ABFF * f AAFF * 0 + f ABFF * f ABFF * 0 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 
f BBFF * f ABEE * 1 / 2 + f BBFF * f BBEE * 1 + f BBFF * f AAEF * 0 + f BBFF * f ABEF * 1 / 4 + 
f BBFF * f BBEF * 1 / 2 + f BBFF * f AAFF * 0 + f BBFF * f ABFF * 0 + f BBFF * f BBFF * 0 + 0] 
f AAFF Drive = FullSimplify [f AAEE * f AAEE * 0 + f AAEE * f ABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 0 + f AAEE * f BBEF * 0 + f AAEE * f AAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * f BBFF * 0 + f ABEE * f AAEE * 0 + f ABEE * f ABEE * 0 + 
f ABEE * f BBEE * 0 + f ABEE * f AAEF * 0 + f ABEE * f ABEF * 0 + f ABEE * f BBEF * 0 + 
f ABEE * f AAFF * 0 + f ABEE * f ABFF * 0 + f ABEE * f BBFF * 0 + f BBEE * f AAEE * 0 + 
f BBEE * f ABEE * 0 + f BBEE * f BBEE * 0 + f BBEE * f AAEF * 0 + f BBEE * f ABEF * 0 + 
f BBEE * f BBEF * 0 + f BBEE * f AAFF * 0 + f BBEE * f ABFF * 0 + f BBEE * f BBFF * 0 + 
f AAEF * f AAEE * 0 + f AAEF * f ABEE * 0 + f AAEF * f BBEE * 0 + f AAEF * f AAEF * 1 / 4 + 
f AAEF * f ABEF * 1 / 8 + f AAEF * f BBEF * 0 + f AAEF * f AAFF * 1 / 2 + 

f AAEF * f ABFF * 1 / 4 + f AAEF * f BBFF * 0 + f ABEF * f AAEE * 0 + f ABEF * f ABEE * 0 + 

f ABEF * f BBEE * 0 + f ABEF * f AAEF * (1 - dl) / 4 + f ABEF * f ABEF * (1 - dl) / 8 + 

f ABEF * f BBEF * 0 + f ABEF * f AAFF * (1 - dl) / 2 + f ABEF * f ABFF * (1 - dl) / 4 + 

f ABEF * f BBFF * 0 + f BBEF * f AAEE * 0 + f BBEF * f ABEE * 0 + f BBEF * f BBEE * 0 + 

f BBEF * f AAEF * 0 + f BBEF * f ABEF * 0 + f BBEF * f BBEF * 0 + f BBEF * f AAFF * 0 + 

f BBEF * f ABFF * 0 + f BBEF * f BBFF * 0 + f AAFF * f AAEE * 0 + f AAFF * f ABEE * 0 + 

f AAFF * f BBEE * 0 + f AAFF * f AAEF * 1 / 2 + f AAFF * f ABEF * 1 / 4 + f AAFF * f BBEF * 0 + 

f AAFF * f AAFF * 1 + f AAFF * f ABFF * 1 / 2 + f AAFF * f BBFF * 0 + f ABFF * f AAEE * 0 + 

f ABFF * fABEE * 0 + f ABFF * f BBEE * 0 + f ABFF * f AAEF * (1 - dl) / 2 + 

f ABFF * f ABEF * (1 - dl) / 4 + f ABFF * f BBEF * 0 + f ABFF * f AAFF * (1 - dl) + 

f ABFF * f ABFF * (1 - dl) / 2 + f ABFF * f BBFF * 0 + f BBFF * f AAEE * 0 + 

f BBFF * fABEE * 0 + f BBFF * f BBEE * 0 + f BBFF * f AAEF * 0 + f BBFF * f ABEF * 0 + 

f BBFF * f BBEF * 0 + f BBFF * f AAFF * 0 + f BBFF * f ABFF * 0 + f BBFF * f BBFF * 0 + 0] 

fABFF Drive = FullSimplify [f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 0 + f AAEE * f BBEF * 0 + f AAEE * f AAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * f BBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * f BBEE * 0 + fABEE * f AAEF * 0 + fABEE * f ABEF * 0 + fABEE * f BBEF * 0 + 
fABEE * f AAFF * 0 + fABEE * f ABFF * 0 + fABEE * f BBFF * 0 + f BBEE * f AAEE * 0 + 
f BBEE * fABEE * 0 + f BBEE * f BBEE * 0 + f BBEE * f AAEF * 0 + f BBEE * f ABEF * 0 + 
f BBEE * f BBEF * 0 + f BBEE * f AAFF * 0 + f BBEE * f ABFF * 0 + f BBEE * f BBFF * 0 + 
f AAEF * f AAEE * 0 + f AAEF * fABEE * 0 + f AAEF * f BBEE * 0 + f AAEF * f AAEF * 0 + 
f AAEF * f ABEF * 1 / 8 + f AAEF * f BBEF * 1 / 4 + f AAEF * f AAFF * 0 + f AAEF * f ABFF * 1 / 4 + 
f AAEF * f BBFF * 1 / 2 + f ABEF * f AAEE * 0 + f ABEF * fABEE * 0 + f ABEF * f BBEE * 0 + 
f ABEF * f AAEF * dl / 4 + f ABEF * f ABEF * 1 / 8 + f ABEF * f BBEF * (1 - dl) / 4 + 
f ABEF * f AAFF * dl / 2 + f ABEF * f ABFF * 1 / 4 + f ABEF * f BBFF * (1 - dl) / 2 + 
f BBEF * f AAEE * 0 + f BBEF * fABEE * 0 + f BBEF * f BBEE * 0 + f BBEF * f AAEF * 1 / 4 + 
f BBEF * f ABEF * 1 / 8 + f BBEF * f BBEF * 0 + f BBEF * f AAFF * 1 / 2 + f BBEF * f ABFF * 1 / 4 + 
f BBEF * f BBFF * 0 + f AAFF * f AAEE * 0 + f AAFF * fABEE * 0 + f AAFF * f BBEE * 0 + 
f AAFF * f AAEF * 0 + f AAFF * f ABEF * 1 / 4 + f AAFF * f BBEF * 1 / 2 + f AAFF * f AAFF * 0 + 
f AAFF * f ABFF * 1 / 2 + f AAFF * f BBFF * 1 + f ABFF * f AAEE * 0 + f ABFF * fABEE * 0 + 
f ABFF * f BBEE * 0 + f ABFF * f AAEF * dl / 2 + f ABFF * f ABEF * 1 / 4 + f ABFF * f BBEF * (1 - dl) / 2 + 
f ABFF * f AAFF * dl + f ABFF * f ABFF * 1 / 2 + f ABFF * f BBFF * (1 - dl) + f BBFF * f AAEE * 0 + 
f BBFF * fABEE * 0 + f BBFF * f BBEE * 0 + f BBFF * f AAEF * 1 / 2 + f BBFF * f ABEF * 1 / 4 + 
f BBFF * f BBEF * 0 + f BBFF * f AAFF * 1 + f BBFF * f ABFF * 1 / 2 + f BBFF * f BBFF * 0 + 0] 

fBBFF Drive = FullSimplify [f AAEE * f AAEE * 0 + f AAEE * fABEE * 0 + f AAEE * f BBEE * 0 + 
f AAEE * f AAEF * 0 + f AAEE * f ABEF * 0 + f AAEE * f BBEF * 0 + f AAEE * f AAFF * 0 + 
f AAEE * f ABFF * 0 + f AAEE * f BBFF * 0 + fABEE * f AAEE * 0 + fABEE * fABEE * 0 + 
fABEE * f BBEE * 0 + fABEE * f AAEF * 0 + fABEE * f ABEF * 0 + fABEE * f BBEF * 0 + 
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1 

- (4 fAAEE + 2 (fAAEF + fABEE) + fABEF) ( 2 fAAEE + fAAEF -(- 1 + dO ) (2 fABEE + fABEF) 
8 

1 

- (fAAEF (2 fABEE + 4 dO fABEE + fABEF + 2 dO fABEF + 8 fBBEE + 4 fBBEF) + 
8 

2 fAAEE ( (2 + 4 dO) fABEE + fABEF + 2 dO fABEF + 8 fBBEE + 4 fBBEF) + 
(2 fABEE + fABEF) ( 2 f ABEE + f ABEF - ( - 3 + 2 dO ) (2 fBBEE + fBBEF) ) ) 



(dO (2 fABEE + fABEF) + 2 fBBEE + fBBEF) (2 fABEE + fABEF + 4 fBBEE + 2 fBBEF) 

(4 fAAEF 2 + 2 fAAEE (4 fAAEF + 8 fAAFF + 3 fABEF - 2 dO fABEF + 2 fABFF + 4 dO fABFF) - 

(2 fABEE + fABEF) ( ( - 6 + 4 dl ) f AAFF + ( - 2 + dO + dl ) ( fABEF + 2 fABFF) ) + 

fAAEF ( 8 fAAFF - 2 ( ( - 3 + 2 dl ) fABEE + (-3 + dO+dl) fABEF + ( - 3 + 2 dO ) fABFF ) ) ) 

(2 fABEE fABEF + fABEF 2 + 4 fABEE fABFF + 2 fABEF fABFF + 3 fABEF fBBEE - 

2 dO fABEF fBBEE + 2 fABFF fBBEE + 4 dO fABFF fBBEE + 3 fABEE fBBEF - 
2 dl fABEE fBBEF + 3 fABEF fBBEF - dO fABEF fBBEF - dl fABEF fBBEF + 3 fABFF fBBEF - 
2 dO fABFF fBBEF + fAAFF ( (2 + 4 dl) fABEE + fABEF + 2 dl fABEF + 8 fBBEE + 4 fBBEF) - 
( - 3 + 2 dl ) (2 fABEE + fABEF ) fBBFF + 

fAAEE (fABEF + 2 dO fABEF + 6 fABFF - 4 dO fABFF + 4 fBBEF + 8 fBBFF) + fAAEF 

(fABEE + 2 dl fABEE + (1 + dO + dl) fABEF + fABFF + 2 dO fABFF + 4 (fBBEE + fBBEF + fBBFF) ) 

(dO (fABEF 2 + 2 fABEE (fABEF + 2 fABFF) + 

4 fABFF (-2 fBBEE + fBBEF) + 2 fABEF (fABFF + 2 fBBEE + fBBEF) ) + 
2 (fABFF (6 fBBEE + fBBEF) + fABEF (fBBEE + fBBEF + fBBFF) + 

(fABEE + 4 fBBEE + 2 fBBEF) (fBBEF + 2 fBBFF) ) + 
dl (2 fABEE + fABEF) (fABEF + 2 (fABFF + fBBEF + 2 fBBFF) ) ) 

(2 fAAEF + 4 fAAFF + fABEF + 2 fABFF) (fAAEF + 2 fAAFF - ( - 1 + dl) (fABEF + 2 fABFF) ) 
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1 

- (fAAEF ( f ABEF + 2 dl f ABEF + 2 f ABFF + 4 dl f ABFF + 4 fBBEF + 8 fBBFF) + 
8 

2 fAAFF ( f ABEF + 2 dl f ABEF + 2 f ABFF + 4 dl f ABFF + 4 fBBEF + 8 fBBFF) + 
(fABEF + 2 fABFF) ( f ABEF + 2 f ABFF - ( - 3 + 2 dl ) ( fBBEF + 2 fBBFF) ) ) 

1 

- (dl (fABEF + 2 fABFF) + fBBEF + 2 fBBFF) (fABEF + 2 (fABFF + fBBEF + 2 fBBFF) ) 
8 

Selection 

wAAEE = wAAEF = wAAFF = 1 ; 
wABEE = wABEF = wABFF = 1 - hs ; 
wBBEE = wBBEF = wBBFF = 1 - s ; 

W = FullSimplify [ (f AAEE Drive wAAEE + f AAEF Drive wAAEF + f AAFF Drive wAAFF 
+ f ABEE Drive wABEE + f ABEF Drive wABEF + f ABFF Drive wABFF 
+ f BBEE Drive wBBEE + f BBEF Drive wBBEF + f BBFF Drive wBBFF 

)]; 

f AAEE S ei = FullSimplify [fAAEE Drive wAAEE / w] ; 
fABEE S ei = FullSimplify [fABEE Drive wABEE / w] ; 
fBBEE Sel = FullSimplify [fBBEE Drive wBBEE / w] ; 
fAAEFgei = FullSimplify [fAAEF Drive wAAEF / w] ; 
fABEFgei = FullSimplify [fABEF Drive wABEF / w] ; 
fBBEFgei = FullSimplify [fBBEF Drive wBBEF / w] ; 
f AAFF Sel = FullSimplify [fAAFF Drive wAAFF / w] ; 
fABFF S ei = FullSimplify [fABFF Drive wABFF / w] ; 
f BBFF Sel = FullSimplify [fBBFF Drive wBBFF / w] ; 

fA Se i = FullSimplify [fAAEE Sel + f AAEF Sel + f AAFF Sel + ( f ABEE Sel + f ABEF Sel + f ABFF Sel ) / 2] 
fB Se i = FullSimplify [fBBEEgei + f BBEF Sel + f BBFF Sel + ( f ABEE Sel + f ABEF Sel + f ABFF Sel ) / 2] 
fE sei = FullSimplify [fAAEE Sel + f ABEE Sel + f BBEE Sel + ( f AAEF Sel + f ABEF Sel + f BBEF Sel ) / 2] 
fFsei = FullSimplify [fAAFFgei + f ABFF Sel + f BBFF Sel + ( f AAEF Sel + f ABEF Sel + f BBEF Sel ) / 2] 
AfA = FullSimplify [fA Sel - fA] ; 
AfB = FullSimplify [fB Se i - fB] ; 
AfE = FullSimplify [fE Sel - fE] ; 
AfF = FullSimplify [fF Se i - fF] ; 

(*AfC=FullSimplify [fCgei-fC] ;*) 

(*Check: do genotype freqs after selection sum to one?*) 
FullSimplify [fA Sel + fB Sel ] 
FullSimplify [fE Se i + fFsei] 
1 



1 
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Analysis 

Invasion of sperm-acting drive modifier unlinked to the driver 

invasionOfUnlinkedModeif ier = 
FullSimplify [w AfF / fF /. randomHaps /. D -> 0 / • sum2one /. minormod] /. 
fF->0(*this assumes no LD and random mating*) 

- (-1 + fB) fB ( (-1 + 2 fB) hs - fB s) e 

(*With some rearangement , 

invasionOfUnlikedModeif ier= -e fA fB(fB s+(fA-fB) hs) ) 
So long as sh<s, [ 

i.e. the driver does not have an underdominant effect on fitness] this is 

guaranteed to be positive when e is negative [i.e. the modifier dampens drive]*) 

Fixation of sperm-acting drive modifier unlinked to the driver 

f ixationOf UnlinkedModeif ier = FullSimplify [ 

FullSimplify [w AfF / (1 - fF) /. randomHaps /. D -> 0 / • sum2one /. minormod] /. fF -> l] 
(*this assumes no LD and random mating*) 

(- 1 + fB) fB (hs - 2 fB hs + fB s) e 

(★Thus the invasion and fixation conditions are equivalent*) 



